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A  numerical  and  experimental  investigation  of  a  system  for  the  micro-cogeneration  of  heat  and  power, 
based  on  a  Stirling  cycle  and  equipped  with  a  flameless  burner,  is  carried  out  with  the  purpose  of  eval¬ 
uating  the  system  performances  with  hydrogen-enriched  fuels.  The  numerical  model  of  the  combustion 
chamber  gave  significant  insight  concerning  the  analysis  and  interpretation  of  the  experimental  mea¬ 
surements;  however  it  required  considerable  efforts,  especially  regarding  the  definition  of  proper  bound¬ 
ary  conditions.  In  particular  a  subroutine  was  developed  in  order  to  couple  the  oxidation  process  to  the 
overall  operation  of  the  micro-cogeneration  unit.  This  procedure  was  proved  to  perform  satisfactory,  pro¬ 
viding  values  of  the  heat  source  for  the  Stirling  cycle  that  lead  to  expected  thermal  efficiencies  in  agree¬ 
ment  with  those  indicated  in  the  literature  for  similar  systems.  The  importance  of  a  proper  turbulence- 
chemistry  interaction  treatment  and  rather  detailed  kinetic  schemes  to  capture  flameless  combustion 
was  also  assessed.  A  simple  NO  formation  mechanism  based  on  the  thermal  and  prompt  routes  was  found 
to  provide  NO  emissions  in  relatively  good  agreement  with  experimental  observations  when  applied  on 
thermo-fluid  dynamic  fields  obtained  from  detailed  oxidation  schemes. 

©  2011  Published  by  Elsevier  Ltd. 


1.  Introduction 

Combined  heat  and  power  generation  (CHP),  or  cogeneration, 
has  recently  received  particular  attention  as  an  alternative  to  con¬ 
ventional  methods  of  producing  heat  and  electricity  separately. 

In  particular,  the  continuous  fluctuations  in  energy  price  and 
the  uncertainty  in  fuel  supplies  have  determined  a  growing  inter¬ 
est  for  the  use  of  micro-cogeneration  in  the  residential  sector.  In 
fact,  micro-CHP  systems  are  characterized  by  very  high  conversion 
efficiencies  (over  80%)  with  respect  to  fossil  fuel  fired  electricity 
generation  systems.  Moreover,  distributed  generation  in  decentral¬ 
ized  systems  is  very  appealing  to  reduce  greenhouse  gases  emis¬ 
sions  and  to  exploit  the  new  opportunities  offered  by  alternative 
fuels,  such  as  hydrogen-based  fuels  obtained  from  biomass  (and 
coal)  thermal  conversion  processes. 

Onovwiona  and  Ugursal  [1]  recently  carried  out  a  literature  re¬ 
view  on  residential  cogeneration  systems,  taking  into  account  four 
different  technologies:  reciprocating  internal  combustion  engines, 
micro-turbines,  fuel  cells  and  Stirling  engines.  The  authors  dis¬ 
cussed  the  development  stage  of  the  technologies,  indicating  envi¬ 
ronmental  benefits,  performances  and  installation  costs.  As  a 
result,  the  privileged  choice  for  residential  application  appeared 
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to  be  the  reciprocating  internal  combustion  engines  (ICE),  due  to 
their  robustness,  reliability  and  reasonable  costs.  However,  the 
“maturity”  of  the  ICE  technology  is  also  responsible  of  the  poor 
performances  of  such  systems  from  the  environmental  point  of 
view,  e.g.  CO  and  NO*  emissions,  although  improvements  have 
been  recently  made  for  micro-CHP  applications  [2].  On  the  other 
hand,  innovative  technologies  such  as  Stirling  engines  are  charac¬ 
terized  by  emissions  about  10  times  smaller  [1]  than  those  of  Otto 
engines  with  catalytic  converters  and  they  offer  more  flexibility 
from  the  point  of  view  of  fuel  nature. 

A  recent  life  cycle  assessment  [3  ]  analysis  on  micro-cogeneration 
has  indicated  micro-CHP  systems  are  superior  in  terms  of  GHG 
emissions  to  both  average  electricity  and  heat  supply  as  well  as  to 
separate  production  of  electricity  in  gas  power  plants  and  heat  in 
condensing  boilers.  Among  the  existing  possibilities  (i.e.  fuel  cells, 
micro-turbines),  cogeneration  systems  based  on  Stirling  cycles  ap¬ 
pear  particularly  appealing,  as  they  ensure  high  efficiency,  fuel  flex¬ 
ibility,  low  emissions,  low  noise/vibration  levels  and  good 
performance  at  partial  load  [1  ].  Since  the  combustion  process  takes 
places  outside  the  engine,  different  fuels  may  be  used;  moreover  the 
combustion  process  can  be  controlled  independently  from  the  en¬ 
gine  operation.  Besides,  the  maintenance  of  Stirling  engines  is  less 
critical  compared  to  other  reciprocating  internal  combustion  en¬ 
gines,  because  of  fewer  moving  parts  leading  to  quieter  and  smooth¬ 
er  operations.  The  versatility  regarding  the  fuel  topology  is 
extremely  important  as  it  implies  that  non  conventional  fuels  such 
as  biomasses  and  low-calorific  gases  enriched  with  hydrogen  could 
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Nomenclature 

Cl 

confidence  interval 

A  Tmax 

maximum  temperature  rise,  I< 

& 

specific  heat,  J  kg-1  K_1 

<'AIR 

air  excess 

E 

error  associated  to  the  numerical  simulations 

S 

kinetic  energy  dissipation  rate,  m2  s-3 

Fs 

safety  factor 

t ]th, St 

thermal  efficiency  of  the  Stirling  engine 

GCI 

Grid  Convergence  Index 

f1  reg 

regenerator  efficiency 

k 

turbulent  kinetic  energy,  m2  s-2 

V 

degrees  of  freedom 

recirculation  degree 

h 

grid  spacing,  m 

Subscripts 

n 

number  of  experimental  replicates 

0 

at  ambient  conditions 

ifi 

mass  flow  rate,  kg  s_1 

air 

air 

PHe 

helium  operating  pressure,  Pa 

b 

base  grid 

Qrem 

heat  removed  in  the  flue  gases/helium  heat  exchanger, 

c 

coarse  grid 

W 

e 

electrical 

Qin 

burner  input  power,  W 

i 

fine  grid 

Qair 

air  preheating  power,  W 

fg 

flue  gases 

Qfg 

power  recoverable  from  the  exhaust  gases,  W 

fg,R 

recirculating  flue  gases 

r 

radial  coordinate,  m 

fuel 

fuel 

1~hc 

grid  coarsening  ratio 

k 

species  index 

rhf 

grid  refining  ratio 

m 

mixed  value,  referred  to  the  mixing  process  of  fresh  and 

S 

sample  standard  deviation 

flue  gases 

T 

temperature,  K 

out 

value  at  chimney 

fa/2, v 

quantile  of  the  measured  population 

th 

thermal 

X 

axial  coordinate,  m 

tot 

overall 

ym 

predicted  variable 

ye 

measured  variable 

Abbreviations 

Y 

mass  fraction 

BC 

boundary  condition 

Ue 

experimental  uncertainty 

CHP 

combined  heat  and  power  generation 

Usver 

solution  uncertainty 

EDC 

Eddy  dissipation  concept 

wel 

electric  power,  W 

HITAC 

high  temperature  air  combustion 

MILD 

moderate  or  intense  low-oxygen  dilution 

Greek  symbols 

PDF 

probability  density  function 

a 

desired  level  of  confidence 

WSGG 

weighted  sum  of  gray  gases 

be  effectively  exploited,  without  damaging  the  engine.  Indeed,  the 
performances  of  the  combustion  system  in  terms  of  efficiency  and 
pollutants  emissions  need  to  be  assessed  when  working  with  non 
conventional  fuels. 

Thomas  [4]  recently  performed  a  benchmark  testing  of  four  mi- 
cro-CHP  units  in  the  range  of  5-9  kW  electric  power,  two  of  them 
based  on  Stirling  engines  and  two  Otto  CHP  units.  The  systems 
were  compared  on  the  basis  on  the  German  environmental  regula¬ 
tion  which  fixes  the  admissible  performances  of  CHP  systems  from 
the  point  of  view  of  electrical,  thermal  and  overall  efficiencies,  NOx 
and  CO  emissions.  Results  of  the  benchmark  clearly  showed  that 
the  only  unit  able  to  meet  all  the  regulation  requirements  was  a 
SOLO  Stirling  unit  [5],  equipped  with  a  flameless  burner. 

The  implementation  of  flameless  burners  in  micro-CHP  units 
represents  a  very  attractive  solution,  due  to  the  great  potentials 
of  the  flameless  technology  to  deal  with  complex  fuels.  Flameless 
combustion  [6],  or  MILD  [7,8]  of  HITAC  [9-11]  has  received  partic¬ 
ular  attention  in  the  last  years  due  to  its  ability  of  combining  high 
combustion  efficiencies  with  extremely  low  pollutant  emissions. 
Such  combustion  regime  is  based  on  the  modification  of  the  tradi¬ 
tional  flame  structure:  the  system  is  driven  towards  homogeneous 
(temperature  and  species)  conditions,  by  means  of  a  massive  recir¬ 
culation  of  exhaust  gases  in  the  reaction  region.  As  a  result,  the  sys¬ 
tem  approaches  perfectly  stirred  reactor  conditions  and 
temperature  increase  due  to  combustion  is  diluted  over  a  wider 
reaction  zone.  The  volumetric  nature  of  the  combustion  process  al¬ 
lows  achieving  very  high  combustion  efficiencies  ensuring  near 
zero  CO  emissions  in  the  exhausts.  Moreover,  NOx  emissions  are  lar¬ 
gely  suppressed,  even  when  large  air  preheating  is  applied,  because 
of  the  reduced  temperature  peaks. 


The  latter  effect  appears  particularly  beneficial  for  controlling 
NOx  formation  and  shows  potentials  for  limiting  the  reactivity  of 
hydrogen-based  fuels.  The  effectiveness  of  flameless  combustion 
with  hydrogen-based  fuels  has  been  recently  demonstrated 
[11-15]. 

The  present  paper  aims  at  investigating  a  flameless  system  for 
the  micro-CHP  of  heat  and  power  based  on  a  Stirling  cycle,  with 
the  purpose  of  evaluating  its  performance  with  hydrogen-enriched 
fuels.  A  joint  experimental  and  modeling  activity  has  been  carried 
out  to  gain  insight  into  the  main  features  of  the  system,  which 
shows  complex  geometry,  due  to  its  industrial  characteristics, 
and  operation,  due  to  the  coupling  between  the  combustion  pro¬ 
cess  and  a  Stirling  cycle.  The  available  experimental  data  are,  then, 
crucial  for  a  constructive  validation  of  the  computational  approach, 
with  particular  attention  to  oxidation  kinetic  schemes  and  pollu¬ 
tant  formation  models. 


2.  The  burner 

The  investigated  micro-CHP  system1  [5],  is  able  of  providing 
2-9  kW  of  electricity  and  22-30  kW  of  thermal  power,  with  an  elec¬ 
trical  efficiency  of  22-24%  and  a  total  efficiency  larger  than  90%.  The 
system  is  equipped  with  a  flameless  burner  for  the  oxidation  of  the 
gaseous  fuel  and  a  finned  heat  exchanger  is  placed  inside  the  com¬ 
bustion  chamber  to  supply  power  to  a  Stirling  cycle  operated  with 
a  different  medium  (i.e.  helium).  The  CHP  unit  is  designed  to  burn 
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natural  gas;  however  the  burner  behavior  with  H2  based  fuels  was 
studied  varying  the  H2  mass  fraction  in  the  fuel  up  to  10%. 

The  system  consists  of  several  sub-assemblies  enclosed  by  an 
overall  housing.  The  cylindrical  and  horizontally  arranged  combus¬ 
tion  chamber  (0.31  m  long)  is  closed  at  the  end  and  delimited  by  an 
external  (steel)  case  with  a  diameter  varying  from  0.021  to 
0.0235  m.  Inside  the  combustion  chamber,  there  is  a  flame  tube 
(Fig.  la)  with  a  diameter  of  0.087  m,  into  which  fresh  air  and  fuel 
are  fed,  and  produce  the  oxidation  process.  The  flue  gases  flow 
along  the  chamber  axis  until  they  reach  the  closed  end  and  reverse 
direction,  where  the  external  enclosure  diameter  is  larger.  Then, 
the  exhausts  cross  the  finned  heat  exchanger,  denoted  as  heater 
(Fig.  la),  which  supplies  heat  to  a  helium  stream  carrying  out  a 
Stirling  cycle.  Finally,  a  part  of  exhaust  gases  is  entrained  into 
the  flame  tube  by  the  incoming  air,  fed  to  the  combustion  chamber 
through  six  nozzles  inclined  by  25°  with  respect  to  the  horizontal 
(Fig.  lb).  The  high  momentum  air  jet  also  entrains  the  fuel,  which 
is  injected  at  low  velocity  through  six  nozzles  inclined  by  65°  with 
respect  to  the  horizontal  (Fig.  lb).  The  recirculation  of  the  flue 
gases  into  the  reaction  region  is  essential  to  achieve  flameless  com¬ 
bustion,  as  it  ensures  the  dilution  of  the  fresh  gases  as  well  as  their 
preheating  above  the  fuel  self  ignition  temperature.  The  fraction  of 
exhausts  which  do  not  recirculate  in  the  reaction  zone  leaves  the 
combustion  chamber,  crossing  a  honeycomb  heat  exchanger  for 
air  preheating.  This  allows  increasing  the  overall  efficiency,  provid¬ 
ing  air  at  about  823  K.  Before  the  chimney,  exhaust  gases  are 
cooled  down  in  a  water-gas  heat  exchanger.  An  additional 


water-water  heat  exchanger  is  also  included  in  the  micro-CHP  sys¬ 
tem,  to  recover  the  heat  released  by  helium.  The  overall  operation 
of  the  micro-CFlP  system  with  the  main  fluid  streams  and  the  heat 
exchanges  is  shown  in  Fig.  2.  The  system  performances  can  be  ad¬ 
justed  by  varying  the  helium  working  pressure  in  the  range  40- 
150  bar,  corresponding  to  an  electrical  output  from  2  to  9  kW. 

Regarding  the  experimental  measurements  available,  five  ther¬ 
mocouples  are  placed  on  the  external  surface  of  the  exhausts/helium 
heat  exchanger  to  measure  the  flue  gas  temperature;  moreover,  the 
temperatures  on  top  of  the  expansion  cylinder  of  the  Stirling  cycle 
and  before  the  chimney  are  also  monitored.  As  for  the  species  com¬ 
position,  the  composition  of  the  exhaust  stream  leaving  the  combus¬ 
tion  chamber  is  characterized  (02,  C02,  CO  and  NO). 

2. 2.  Stirling  engine 

As  mentioned  previously,  the  finned  heat  exchanger  (Fig.  1) 
provides  heat  to  helium,  which  drives  a  Stirling  engine  (Fig.  2) 
for  the  generation  of  electrical  power.  The  Stirling  engine  consists 
of  a  compression  and  an  expansion  cylinders,  placed  at  90°.  Heater, 
regenerator  and  gas  cooler  are  positioned  between  the  two  cylin¬ 
ders.  The  working  gas  is  moved  between  the  cylinders,  undergoing 
a  closed  thermodynamic  cycle.  Inside  the  compression  cylinder  the 
gas  is  isothermally  compressed  at  low  temperature  by  water  cool¬ 
ing;  then,  it  is  moved  through  a  regenerator  towards  the  expansion 
cylinder.  The  regenerator,  hot  in  this  case,  heats  the  helium  stream 
at  temperature  levels  of  about  923  K.  During  the  isothermal 


Fig.  1.  Flame  tube  and  finned  heat  exchanger  (a);  air  and  fuel  feeding  system  (b). 


Compression 

cylinder 


Fig.  2.  Schematic  diagram  of  flow  streams  and  heat  exchanges  with  highlighted  combustion  section. 
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expansion  the  gas  is  constantly  heated  by  the  heater;  afterwards, 
the  gas  is  moved  back  through  the  regenerator  (cold  in  this  case), 
where  it  is  cooled  down,  to  the  compression  cylinder.  The  use  of 
a  regenerator  considerably  improves  the  cycle  efficiency. 

3.  Experimental  campaign 

The  experimental  campaign  of  the  micro-CHP  Stirling  unit  has 
been  carried  with  different  operating  conditions,  varying  the 
hydrogen  content  in  the  fuel  stream  and  the  operating  pressure 
of  helium.  Table  1  reports  the  boundary  conditions  for  the  experi¬ 
mental  runs,  whereas  the  main  experimental  results  are  listed  in 
Table  2.  Columns  2-5  in  Table  2  list  the  composition  of  the  flue 
gases  for  the  different  investigated  runs.  It  can  be  observed  that 
oxygen  content  in  the  flue  gases  increases  with  the  hydrogen  mass 
fraction.  This  is  expected  as  the  normal  operation  of  the  system 
leads  to  larger  air  excess  for  hydrogen-enriched  mixtures,  as  indi¬ 
cated  in  Table  1.  In  addition  CO  and  C02  levels  decrease  with  H2, 
because  of  the  lower  relative  carbon  fraction  in  the  fuel.  Finally, 
regarding  NO  emissions,  a  slight  decrease  is  observed  with  increas¬ 
ing  H2,  for  higher  burner  loads  (i.e.  helium  pressures,  PHe,  of  90  and 
120  bar).  This  is  due  to  the  increase  of  air  excess  when  hydrogen  is 
added  to  the  fuel,  which  leads  to  a  higher  dilution  of  the  flue  gases. 


Table  1 

Summary  of  the  experimental  campaign  on  the  the  CHP  unit. 


Run 

mNG  (kg/s) 

mH2  (kg/s) 

Yh2  (-) 

£air  (%) 

Qi„  (1<W) 

PHe 

1 

5.23E-04 

0 

0 

48 

26.1 

90 

2 

4.47E-04 

3.14E-05 

0.07 

73 

26.2 

90 

3 

4.42E-04 

4.34E-05 

0.09 

82 

27.3 

90 

4 

6.39E-04 

0 

0 

43 

32.0 

120 

5 

5.62E-04 

3.97E-05 

0.07 

69 

32.9 

120 

6 

5.37E-04 

5.62E-05 

0.09 

79 

33.6 

120 

Columns  6-10  of  Table  2  summarize  the  performances  of  the 
system  at  the  investigated  conditions.  The  electrical,  Wd,  and  ther¬ 
mal  power,  Qffr,  are  listed,  together  with  the  efficiencies  (electric, 
rje,  thermal,  rjthl  and  overall,  ^fot).  It  can  be  observed  how  all  the 
efficiencies  rise  with  the  burner  load,  i.e.  helium  pressure,  reaching 
values  around  97-98%.  At  lower  burner  loads  the  system  perfor¬ 
mances  are  lower,  as  expected,  because  the  system  operates  in 
off-design  conditions  and  the  influence  of  the  heat  losses  is  more 
significant. 

4.  Numerical  modeling  of  the  micro-CHP  flameless  unit 

The  numerical  modeling  of  the  Stirling  CHP  unit  has  been  car¬ 
ried  out  with  the  commercial  software  FLUENT  6.3  by  Ansys  Inc. 

The  activity  has  been  focused  on  the  combustion  section  of  the 
system.  Hypotheses  regarding  the  overall  operation  of  the  micro- 
CHP  unit  have  been  made  to  determine  the  heat  exchange  between 
the  flue  gases  and  the  operating  fluid,  i.e.  helium,  of  the  Stirling  cy¬ 
cle.  The  following  sections  describe  the  approach  adopted,  which  is 
also  summarized  in  the  scheme  of  Fig.  3. 

4.1.  Computational  domain  and  grid 

The  geometry  of  the  combustion  chamber  is  rather  complex, 
due  to  the  presence  of  six  injection  nozzles  for  both  air  and  fuel 
streams  (Fig.  4).  Moreover,  a  detailed  description  of  the  flue 
gases/helium  heat  exchanger  (Fig.  la)  is  not  feasible,  due  to  the 
large  number  of  grid  cells  that  would  be  required  to  model  the  geo¬ 
metrical  details  of  the  heater.  Simplifications  need  to  be  made,  to 
minimize  the  impact  of  the  mesh  on  the  overall  computational  re¬ 
sources,  while  ensuring  an  acceptable  accuracy  in  the  representa¬ 
tion.  Following  the  above  considerations,  it  has  been  chosen  to 
model  a  60°  angular  section  of  the  real  combustion  chamber.  The 


Table  2 

Flue  gas  composition  (dry  basis)  and  main  operation  outputs  for  the  runs  investigated  in  Table  1. 


Run 

CO  (ppmv) 

NO  (ppmv) 

02  (%  vol) 

C02  (%  vol) 

Wel  (1<W) 

Qth  (kW) 

4e(~) 

4th  (-) 

4  tot  (  ) 

1 

12.0 

50.5 

6.4 

8.1 

7.0 

17.0 

0.65 

0.27 

0.92 

2 

6.2 

45.6 

8.5 

6.1 

7.2 

17.6 

0.67 

0.28 

0.95 

3 

5.9 

45.8 

7.8 

6.7 

7.3 

17.8 

0.65 

0.27 

0.92 

4 

17.4 

49.0 

5.9 

8.3 

9.1 

22.0 

0.69 

0.28 

0.97 

5 

7.9 

40.8 

8.0 

6.4 

9.1 

23.0 

0.70 

0.28 

0.98 

6 

6.6 

39.0 

8.6 

5.8 

9.1 

23.5 

0.70 

0.27 

0.97 

Subroutines 


CFD  model  of  the 
combustion  chamber 


Air  Pre-heater 
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Fig.  3.  Scheme  of  the  modeling  approach. 
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Fig.  4.  Sketch  of  the  combustion  chamber  and  details  about  the  feeding  system. 


heat  exchanger  has  been  treated  as  a  fluid  sub-domain  with  a  heat 
sink  equal  to  the  heat  transferred  from  the  flue  gases  to  helium. 
Details  about  the  calculation  of  the  heat  removed  by  the  heater 
are  given  in  Section  4.2.  Finally,  the  flame  tube  (Fig.  la)  has  been 
modeled  as  a  solid  domain  placed  inside  the  combustion  chamber. 

A  hybrid  approach  has  been  adopted  to  reduce  the  overall  num¬ 
ber  of  computational  cells  (Fig.  5).  For  the  fluid  domain,  a  hexahe- 
dral  mesh  has  been  used  everywhere  with  the  exception  for  the 
injection  region,  where  a  tetrahedral  mesh  is  preferred  to  repro¬ 
duce  accurately  injection  nozzles.  The  air  and  injection  nozzles 
have  been  extended  upstream  the  nozzles  of  about  five  character¬ 
istic  diameters,  to  ensure  the  flow  to  be  fully  developed  at  the  bur¬ 
ner  exit  plane.  The  cell  size  is  not  uniform  and  a  higher  density  is 


observable  in  the  injection  region,  to  provide  a  better  characteriza¬ 
tion  of  the  reaction  structure.  The  flame  tube  is  also  meshed  with 
hexahedra.  The  resulting  total  number  of  cells  for  the  base  mesh 
is  about  382.000.  A  grid  sensitivity  analysis  has  been  carried  out 
to  estimate  the  solution  error  associated  to  the  discretization 
adopted.  To  this  purpose,  two  additional  grids  have  been  taken  into 
account,  a  finer  and  a  coarser  one.  The  finer  grid  has  been  obtained 
by  applying  a  uniform  refining  ratio,  i.e.  rhf=  hf\hh  =  0.8,  to  the  base 
mesh,  in  axial  and  radial  directions  Similarly,  the  coarser  grid  is  ob¬ 
tained  by  applying  a  constant  coarsening  factor  rhc  =  \ic\hh  =  1.25  to 
the  base  mesh.  The  total  number  of  grid  cells  for  the  finer  and  coar¬ 
ser  grid  is  550.000  and  270.000,  respectively.  The  meshes  have 
been  compared  in  a  cold  flow  simulation,  with  boundary 


Fig.  5.  Computational  domain  and  grid  with  details  of  the  injection  nozzles. 
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conditions  taken  from  run  4  in  Table  2.  A  solution  uncertainty, 
Usven  has  been  estimated  following  Logan  and  Nitta  approach 
[16].  The  Grid  Convergence  Index  (GCJ)  has  been  also  estimated, 
assuming  a  safety  factor,  Fs  =  1.25,  as  advised  by  Roache  [17]  for 
a  set  of  three  meshes.  Table  3  lists  the  values  of  Usver  and  GCI  for 
the  set  of  investigated  meshes.  The  values  have  been  obtained  by 
averaging  the  radial  profiles  of  axial  velocity  within  the  flame  tube. 
Table  3  clearly  indicates  values  of  Usver  and  GCI  for  the  base  grid 
below  5-7%,  which  can  be  considered  a  reasonable  compromise 
between  CPU  time  requirements  and  the  numerical  accuracy. 

4.2.  Boundary  conditions 

As  far  as  boundary  conditions  (BCs)  are  concerned,  mass  flow  rate 
BCs  are  used  to  specify  the  inlets.  Fuels,  i.e.  methane  and  hydrogen, 
are  fed  at  room  temperature  whereas  air  temperature  downstream 
of  the  pre-heater  is  not  known  from  the  experimental  measure¬ 
ments.  Therefore,  it  has  been  evaluated  from  the  exhaust  gas  tem¬ 
perature  by  assuming  a  preheating  efficiency  of  50%.  Such  value 
ensures  a  preheating  temperature  of  about  823  K  for  nominal  oper¬ 
ating  conditions,  with  the  burner  fed  with  natural  gas  and  a  helium 
pressure  of  120  bar,  as  indicated  by  the  burner  manufacturer  [5]. 

An  user-defined  subroutine  has  been  developed  to  model  the 
heat  exchanger  behavior  and  to  estimate  the  air  inlet  temperature 
(Fig.  3).  The  heater  is  modeled  as  a  fluid  domain  with  a  heat  sink. 
Therefore,  the  external  subroutine  performs  an  energy  balance 
on  the  overall  unit  to  determine  the  heat  removed  by  the  exhaust 
stream.  In  particular,  the  following  operations  are  performed: 

1.  The  average  temperature  of  the  exhaust  gases  is  calculated,  by 
averaging  the  temperatures  of  each  cell  adjacent  to  the  outlet 
boundary  condition. 

2.  Assuming  a  preheating  efficiency  equal  to  50%,  the  air  inlet  tem¬ 
perature  is  estimated  from  the  exhaust  gas  temperature. 

3.  The  heat  recoverable  from  the  exhaust  gases,  Qfg,  and  the  air 
preheating  heat,  Qa!r,  are  calculated,  given  the  flue  gas  temper¬ 
ature  at  the  chimney,  7}&0lit,  and  the  ambient  air  temperature, 
Taij-o,  known  from  the  experimental  campaign  (Table  1): 

Qfg  =  ™fg cp fg(Tfg)(Tfg  -  Tfg  out)  ^  x 

Qair  =  mairCp,air(Tair)(Tair  ~  T'airp) 

where  cpfg(Tfg)  and  Cp,air  ( Tair )  are  the  specific  heat  at  constant 
pressure  calculated  at  average  flue  and  air  temperatures,  respec¬ 
tively.  The  coefficients  for  the  evaluation  of  the  specific  heat  are 
taken  from  the  CHEMKIN  Thermodynamic  Database  [18]. 

4.  An  energy  balance  is  carried  out  on  the  control  region  indicated 
in  Fig.  2,  to  determine  the  heat  transferred  from  the  flue  gases 
to  helium  in  the  combustion  chamber.  Then,  the  heat  removed 
in  the  flue  gases/helium  heat  exchanger,  Qrem,  is: 

Qrem  —  Qin  ~  Qfg  "T  Qair  (2) 

where  Qin  is  the  thermal  input,  known  from  the  experimental 
campaign  (Table  1). 

5.  The  heat,  Qrem,  calculated  at  step  4  is  specified  as  a  volumetric 
heat  sink  for  the  fluid  region  corresponding  to  the  flue  gases/ 
helium  heat  exchanger  (Fig.  4),  thus  cooling  down  the  flue  gases 
crossing  such  fluid  region  accordingly, 


Table  3 

Solution  uncertainty  and  Grid  Convergence  Index  for  the  fine,  base  and  coarse 
computational  grids. 


Mesh 

Elements 

Usver  {%) 

GCI  (%) 

Coarse 

270,000 

7 

9 

Base 

382,000 

5 

7 

Fine 

550,000 

4 

5 

Regarding  the  other  boundary  conditions,  standard  wall  func¬ 
tions  are  used  at  the  walls  and  symmetry  conditions  are  applied 
to  the  planes  delimiting  the  modeled  angular  sector. 

4.3.  Physical  models 

Favre-averaged  Navier-Stokes  equations  are  solved  across  the 
computational  domain,  using  the  standard  k-s  turbulence  model 
for  Reynolds  stresses.  Following  the  results  obtained  in  previous 
studies  on  the  mathematical  modeling  of  flameless  combustion 
[19-21],  the  Eddy  Dissipation  Concept  (EDC)  [22,23]  model  is  ta¬ 
ken  into  account  for  turbulence  chemistry-interactions,  to  capture 
the  volumetric  features  of  the  investigated  combustion  regime.  The 
effect  of  the  kinetic  mechanism  on  the  results  is  investigated,  com¬ 
paring  the  results  provided  by  a  global  kinetic  approach  and  a  de¬ 
tailed  kinetic  scheme.  The  4-step  global  kinetic  scheme  by  Jones 
and  Lindstedt  [24]  is  used  for  both  methane  and  hydrogen-en¬ 
riched  fuels.  This  scheme  will  be  denoted  as  JL.  The  detailed 
scheme  is  the  so-called  KEE-58  [25],  which  involves  17  species  in 
58  reversible  reactions. 

As  far  as  NO  formation  is  concerned,  only  thermal  and  prompt 
NO  formation  are  considered,  due  to  the  relatively  high  process 
temperatures  observed  experimentally,  which  make  the  contribu¬ 
tions  of  other  routes,  i.e.  N20  [26]  and  NNH  [27-29],  negligible,  as 
shown  in  [21].  Thermal  NO  formation  is  modeled  with  a  simplified 
one-step  kinetic  mechanism,  available  in  the  code  [30]  and  ob¬ 
tained  from  the  reduction  of  the  full  Zeldovich  scheme  by  assum¬ 
ing  a  steady  state  for  the  N  radicals  and  relating  the  O  radical 
concentration  to  that  of  oxygen  by  means  of  the  dissociation  reac¬ 
tion.  The  prompt  NO  formation  from  methane  is  modeled  using  a 
similar  approach,  according  to  the  one-step  mechanism  proposed 
by  [31].  For  both  thermal  and  prompt  NO  formations,  the  Arrhe¬ 
nius  equation  is  integrated  over  a  probability  density  function 
(PDF)  for  temperature,  to  account  for  turbulence  effects. 

Finally,  radiation  modeled  with  the  Pi  radiation  model,  based 
on  the  expansion  of  the  radiation  intensity  into  an  orthogonal  ser¬ 
ies  of  spherical  harmonics  [32].  The  Pi  model  can  provide  reason¬ 
able  results  when  radiation  is  emitted  isotropically,  as  it  happens 
for  combustion  applications.  The  Weighted-Sum-of-Gray-Gases 
(WSGG)  model  is  used  for  predicting  the  spectral  properties  [33]. 

The  influence  of  molecular  diffusion  has  been  also  considered. 
Binary  diffusion  coefficients  are  calculated  following  the  kinetic 
theory  and  a  modification  of  the  Chapman-Enskog  formula  [34]. 
Then,  an  effective  diffusion  coefficient  of  the  species  in  the  mixture 
is  obtained  by  applying  the  Wilke’s  mixing  rule  [35]. 

5.  Results 

The  results  of  the  numerical  modeling  on  the  CHP  flameless  are 
shown  here.  First,  the  main  features  of  the  flow  field  inside  the 
combustion  chamber  are  presented.  Then,  the  effect  of  the  kinetic 
mechanism  adopted  for  gas-phase  oxidation  is  analyzed,  followed 
by  a  discussion  regarding  the  influence  of  hydrogen  addition  to 
the  structure  of  the  reaction  zone  in  the  combustion  chamber.  Fi¬ 
nally,  the  validation  of  the  computational  approach  is  shown, 
comparing  the  numerical  results  with  the  available  experimental 
data. 

5.1.  Flow  field  and  main  operating  features 

The  analysis  of  the  flow  field  in  the  combustion  chamber  indi¬ 
cates  that  the  flue  gases  are  entrained  into  the  flame  tube  by  the 
high  velocity  air  jet.  Fuel  is  also  entrained  by  air  and  mixed  with 
flue  gas  stream.  The  recirculation  degree  in  the  combustion  cham¬ 
ber  can  be  determined  as: 
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(til fuel  “1“  Wlair) 


where  mfg.R  is  the  mass  flow  rate  of  the  exhaust  gases  recirculating 
into  the  reaction  zone,  whereas  mfuei  and  mair  represent  the  fuel  and 
air  mass  flow  rates  fed  to  the  burner.  The  values  of  kR  obtained  for 
the  different  runs  are  listed  in  the  second  column  of  Table  4:  the  bur¬ 
ner  generates  a  degree  of  recirculation,  i.e.  106-112%.  The  average 
conditions  in  terms  of  oxygen  availability  (Column  3,  Table  4),  and 
temperature,  Tm  (Column  4,  Table  4),  after  the  flue  and  fresh  gas  mix¬ 
ing  process,  are  approximately  Y02,m  =  0.15  and  Tm  =  1240  K,  respec¬ 
tively.  Therefore,  recirculation  allows  achieving  both  dilution  and 
preheating  of  the  reactants,  required  to  operate  in  flameless  regime. 
The  fifth  column  in  Table  4  lists  the  maximum  temperature  increase 
observed  for  the  different  investigated  runs,  as  provided  by  the  KEE- 
58  kinetic  mechanism.  Following  the  criterion  of  Cavaliere  and  de 
Joannon  [7],  the  burner  operates  in  flameless  mode  only  when  fed 
with  CH4,  being  the  observed  A Tmax  below  the  fuel  self-ignition  tem¬ 
perature  (i.e.  853  K).  When  hydrogen  is  added  to  the  fuel,  the  self¬ 
ignition  temperature  is  in  between  that  of  hydrogen  (i.e.  783  K)  and 
that  of  methane  (i.e.  853  K).  In  such  cases  the  observed  ATmax  values 
are  higher,  i.e.  ~1 020-1 050  K,  thus  indicating  a  strong  effect  of 
hydrogen  reactivity  on  the  flame  structure.  Therefore,  the  system 
operation  cannot  be  strictly  regarded  as  flameless  and  higher  recircu¬ 
lation  degrees  should  be  ensured  when  operating  with  hydrogen- 


Table  4 

Macro  indicators  of  the  main  operating  features  of  the  micro-CHP  unit.  Data  extracted 
from  the  numerical  simulations  carried  out  with  the  KEE-58  kinetic  mechanism. 


Run 

1<R  (%) 

Vm  (-) 

Tm  (I<) 

A Tmax  (I<) 

W el  /  Q.rem  (  ) 

1 

107 

0.15 

1231 

817 

0.45 

2 

110 

0.15 

1234 

1057 

0.48 

3 

112 

0.16 

1251 

1059 

0.47 

4 

106 

0.14 

1196 

675 

0.44 

5 

108 

0.16 

1203 

1016 

0.47 

6 

109 

0.16 

1205 

1022 

0.46 

based  fuels.  The  last  column  in  Table  4  reports  the  ratio  Wei/Qrem  ob¬ 
tained  from  the  numerical  simulations,  using  the  user-defined  func¬ 
tion  described  in  Section  4.2.  Such  ratio  is  related  to  the  thermal 
efficiency  of  the  Stirling  cycle  by  the  regenerator  efficiency: 


^ 1th  ,St 


Wei 

Qrem  +  (1  hreg)Qreg 


where  Qrem  and  r/reg  are  the  heat  supplied  to  the  working  fluid  in  the 
regenerator  and  the  regenerator  efficiency,  respectively.  In  real  life 
operation,  r/reg  is  smaller  than  1  and,  therefore,  the  thermal  Stirling 
efficiency,  rjthjSt  will  be  lower  than  the  ratio  Wei/Qrem .  The  experi¬ 
mental  measurements  collected  in  the  present  investigation  do 
not  allow  the  calculation  of  Qrem  And  the  estimation  of  r\reg.  How¬ 
ever,  the  measured  Wei/Qrem  ratios  (column  6  in  Table  4),  in  the 
range  0.44-0.48,  result  compatible  with  the  expected  performances 
of  a  Stirling  engine  with  regenerator,  which  are  reported  to  be  in  the 
range  35-40%.  Such  result  can  be  considered  as  a  first  validation,  of 
the  model  adopted  to  replace  the  explicit  operation  of  the  finned 
helium/flue  gases  heat  exchanger. 


5.2.  Effect  of  the  kinetic  mechanism  on  the  flame  structure 

Fig.  6  shows  the  contour  plots  of  temperature  inside  the  com¬ 
bustion  chamber  fed  with  CH4,  at  a  helium  pressure,  PHe  =  90  bar, 
predicted  by  the  EDC  model  with  the  JL  (a)  and  KEE-58  (b)  kinetic 
mechanisms.  A  major  influence  of  the  kinetic  scheme  on  the  struc¬ 
ture  of  the  reaction  zone  is  observed.  The  global  chemistry  ap¬ 
proach  results  in  a  flame  shape  characteristic  of  a  non-premixed 
flame  with  a  thin  flame  front  and  high  temperature  peaks  along 
the  burner  axis.  Conversely,  when  the  KEE-58  scheme  is  employed, 
the  reaction  region  is  fully  detached  from  the  burner  nozzles  and 
shifted  along  the  axis,  towards  the  closed  end  of  the  combustion 
chamber.  Moreover  the  temperature  distribution  in  the  flame  tube 
is  much  more  uniform  and  the  maximum  predicted  temperature, 
i.e.  2050  K,  is  significantly  smaller  than  that  provided  by  the  4-step 
(JL)  mechanism,  i.e.  2585  K. 


Fig.  6.  Temperature  distribution  in  the  combustion  chamber  fed  with  CH4,  predicted  by  the  JL  (a)  and  KEE-58  (b)  kinetic  mechanisms.  PHe  =  90  bar.  Run  1,  Table  1. 
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Fig.  7.  Temperature  distribution  in  the  combustion  chamber  fed  with  CH4/H2  (H2  9%  by  wt.)  predicted  by  the  JL  (a)  and  KEE-58  (b)  kinetic  mechanisms.  PHe  =  90  bar.  Run  1, 
Table  1. 


Fig.  7  illustrates  the  temperature  distribution  in  the  combustion 
chamber  fed  with  a  CH4/H2  mixture  containing  9%  (by  wt.)  of  H2,  at 
a  helium  operating  pressure  PHe  =  90  bar.  When  hydrogen  is  added 
to  the  fuel,  the  qualitative  differences  between  the  results  provided 
by  the  JL  and  KEE-58  kinetic  mechanism  are  smaller.  Both  schemes 
indicate  the  ignition  of  the  fuel  mixture  in  the  region  immediately 
downstream  the  fuel  injection  nozzles,  with  a  defined  flame  region 
developed  along  the  burner  axis.  However,  the  results  provided  by 
KEE-58  indicate  a  more  diffuse  reaction  zone,  without  the  thin 


flame  sheets  characterizing  the  temperature  distribution  predicted 
by  the  global  approach.  Finally,  the  temperature  levels  are  signifi¬ 
cantly  lower  with  KEE-58,  being  the  maximum  predicted  temper¬ 
ature,  i.e.  2309  K,  approximately  350  K  lower  than  that  obtained 
with  the  JL  mechanism. 

Such  considerations  are  quantitatively  supported  by  the  analysis 
of  the  radial  profile  of  temperature  provided  by  the  JL  and  KEE-58 
mechanisms,  shown  in  Fig.  8  for  different  axial  locations,  i.e. 
x  =  0.075  m  (a),  x  =  0.150  m  (b),  x  =  0.250  m  (c)  and  x  =  0.350  m  (d). 


Fig.  8.  Radial  temperature  profiles  at  different  axial  locations,  i.e.  x  =  0.075  m  (a),  x  =  0.150  m  (b),  x  =  0.250  m  (c)  and  x  =  0.350  m  (d),  predicted  by  the  JL  (JL)  and  KEE-58 
mechanisms  for  a  hydrogen  content  in  the  fuel  equal  to  7%  and  9%  (by  wt.),  respectively.  PHe  =  90  bar.  Runs  2  and  3,  Table  1. 
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Two  cases,  corresponding  to  7%  and  9%  H2  mass  fractions  are  consid¬ 
ered,  at  a  helium  pressure  PHe  =  90  bar.  The  profiles  provided  by  the 
global  chemistry  approach  not  only  indicate  higher  temperature  lev¬ 
els  with  respect  to  KEE-58,  but  they  also  show  a  different  flame 
structure.  With  KEE-58,  the  temperature  distribution  in  the  high 
temperature  region  (x  =  0.075  m  and  x  =  0.150  m)  is  always  flatter 
than  that  obtained  with  the  JL  scheme,  which  provides  more  pro¬ 
nounced  peaks  indicating  the  existence  of  a  thin  diffusive  flame 
front  (Fig.  8a  and  b).  When  the  distance  from  the  burner  is  increased 
(x  =  0.250  m  and  x  =  0.350  m),  the  profiles  obtained  with  the  two  ap¬ 
proaches  are  qualitatively  very  similar,  although  global  chemistry 
results  in  higher  temperatures. 

5.3.  Effect  of  H2  addition  to  the  fuel 

Fig.  9  shows  the  effect  of  hydrogen  addition,  i.e.  0  (a),  7  (b)  and  9 
(c)  %  (by  wt.),  on  the  flame  structure  in  the  combustion  chamber,  at 
a  helium  pressure  ^He  =  90  bar.  Hydrogen  has  a  significant  impact 
on  the  radical  pool  and,  then,  it  directly  affects  the  ignition  process. 
For  the  methane  case  (Fig.  9a),  the  oxidation  occurs  close  to  the 
exit  section  of  the  flame  tube,  with  a  homogeneous  temperature 
distribution.  When  hydrogen  is  added  to  the  fuel  (Fig.  9b  and  c), 
the  ignition  zone  is  shifted  backward  and  a  well-defined  high  tem¬ 
perature  region  is  observed.  Moreover,  the  flame  becomes  shorter 
and  thinner,  when  increasing  the  H2  mass  fractions.  The  visual  evi¬ 


dence  provided  by  Fig.  9  supports  the  conclusions  suggested  by  the 
ATmax  values,  discussed  in  Section  5.1.  When  H2  is  added  to  the 
fuel,  a  homogeneous  temperature  distribution  is  no  longer  ob¬ 
served  in  the  combustion  chamber,  being  the  effect  of  hydrogen 
reactivity  dominant  with  respect  to  flue  gas  dilution.  A  possible 
solution  to  such  problem  could  imply  the  modification  of  the  bur¬ 
ner  feeding  system,  to  attain  higher  values  of  kR  within  the  com¬ 
bustion  chamber. 

5.4.  Effect  of  burner  load 

Fig.  10  shows  the  radial  temperature  profiles  obtained  at  differ¬ 
ent  axial  locations,  i.e.  x  =  0.075  m  (a),  x  =  0.150  m  (b),  x  =  0.250  m 
(c)  and  x  =  0.350  m  (d),  when  increasing  the  H2  mass  fraction  in  the 
fuel,  for  two  different  helium  pressures,  PHe  =  90  and  120  bar,  cor¬ 
responding  to  burner  loads  of  ~24  and  ~30  kW,  respectively. 
Fig.  10  quantitatively  confirms  the  effect  of  hydrogen  on  fuel  igni¬ 
tion,  as  discussed  in  Section  5.4.  Regarding  the  comparison  be¬ 
tween  the  different  operating  pressures,  it  could  appear 
unexpected  that  the  lower  burner  load  results  in  higher  tempera¬ 
ture  levels,  especially  closer  to  the  burner  (Fig.  10a  and  b).  How¬ 
ever,  such  behavior  can  be  easily  explained  by  taking  into 
account  that  the  burner  is  designed  for  the  highest  load,  so  that 
the  internal  fluid  dynamic  is  optimized  for  such  conditions;  when 
operating  with  reduced  flow  rates,  ignition  occurs  before  the  fresh 


Fig.  9.  Temperature  distribution  in  the  combustion  chamber  fed  with  mixtures  containing  0  (a),  7  (b)  and  9  (c)  %  (by  wt.)  of  H2,  predicted  by  the  KEE-58  kinetic  mechanism. 
PHe  =  90  bar.  Runs  2  and  3,  Table  1. 
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Fig.  10.  Radial  temperature  profiles  at  different  axial  locations,  i.e.  x  =  0.075  m  (a),  x  =  0.150  m  (b),  x  =  0.250  m  (c)  and  x  =  0.350  m  (d),  predicted  by  the  KEE-58  kinetic 
mechanism  for  a  H2  content  in  the  fuel  equal  to  0%,  7%  and  9%  (by  wt.),  respectively.  PHe  =  90  and  120  bar.  Runs  1-6,  Table  1. 


gases  have  entrained  a  sufficient  amount  of  flue  gases,  thus  result¬ 
ing  in  higher  temperature  peaks.  Fig.  10  also  points  out  that  the 
reaction  zone  is  shifted  towards  the  end  of  the  flame  tube  for  the 
CH4  fuel  case.  For  PHe  =  90  bar,  the  reaction  region  is  collocated  be¬ 
tween  x  =  0.250  m  (Fig.  10c)  andx  =  0.350  m  (Fig.  lOd),  as  indicated 
by  the  temperature  rise  at  the  axis,  from  about  1250  I<  to  approx¬ 
imately  2050  K.  On  the  other  hand,  for  the  PHe  =120  bar  case,  no 
ignition  is  observed  up  to  an  axial  distance  of  x  =  0.350  m 
(Fig.  lOd),  the  temperature  profile  simply  increasing  from  the 
mixed  temperature  value,  Tm  (Column  4,  Table  4)  to  the  flame  tube 
temperature.  Then,  the  reaction  zone  is  shifted  further  towards  the 
closed  end  of  the  combustion  chamber  with  respect  to  the 
PHe  =  90  bar  case. 


5.5.  Model  validation 


A  first  validation  of  the  modeling  approach  has  been  already 
provided  in  Section  5.1,  being  the  estimated  heat  flux  to  the  work¬ 
ing  fluid  of  the  Stirling  engine  compatible  with  the  values  of  ther¬ 
mal  efficiencies  found  in  the  literature  for  similar  systems  [36]. 

The  present  section  provides  a  comparison  between  predicted 
and  measured  CO  and  NO  emissions  from  the  unit.  Following 
Oberkampf  and  Barone  [37]  and  Parente  et  al.  [38],  the  error  asso¬ 
ciated  to  the  numerical  simulations,  E,  and  the  related  confidence, 
Cl,  are  determined: 

E  =  (ym-ye)±CI  (5) 


Cl 


where  ym  and  ye  represent  the  model  output  and  the  measured  va¬ 
lue,  respectively  and  n  is  the  number  of  experimental  replicates.  The 
terms  under  the  square  root  in  Eq.  (6)  represent  the  two  contribu¬ 
tions  to  the  global  uncertainty,  the  statistical  reproducibility  and 
the  experimental  error,  respectively.  The  first  term  is  calculated 
assuming  that  the  measurements  follow  a  t-student  distribution: 
ta/ 2>v  is  the  quantile  of  the  measured  population  with  v  =  n  -  1  de¬ 


grees  of  freedom,  s  is  the  sample  standard  deviation  and  a  is  the  de¬ 
sired  level  of  confidence.  The  experimental  uncertainty,  Ue,  is 
obtained  from  the  literature.  Following  Gluck  et  al.  [39]  and  Moun- 
aud  [40],  relative  errors  of  5%  and  3%  have  been  assumed  for  the 
measurements  of  NO  and  CO,  carried  out  with  the  Siemens  ULTRA¬ 
MAT  6  analyzer. 

Global  metrics  can  be  also  provided  to  summarize  the  perfor¬ 
mances  of  the  two  kinetic  mechanisms  for  the  two  investigated 
pressure  levels,  PHe  =  90  and  120  bar.  In  particular,  an  average  error 
metric  and  confidence  interval  can  be  defined: 


Fig.  11  shows  the  comparison  between  measured  and  calcu¬ 
lated  CO  (a)  and  NO  (b)  emissions  from  the  combustion  chamber, 
with  the  different  kinetic  mechanisms  adopted.  It  is  very  clear  that 
the  global  kinetic  approach  results  in  major  discrepancies  for  both 
CO  and  NO  values,  indicating  its  inadequacy  to  properly  describe 
the  oxidation  process;  in  particular,  the  large  NO  emissions  reflect 
a  non  negligible  over  prediction  of  the  temperature  levels  by  the  JL 
scheme,  as  already  discussed  in  Section  5.3.  With  regard  to  the 
KEE-58  kinetic  mechanism,  a  very  good  agreement  with  the  exper¬ 
imental  data  is  observed  for  the  =  90  bar  case.  An  excellent 
quantitative  agreement  is  observed  for  both  NO  and  CO  emissions, 
as  confirmed  by  Table  5.  The  average  error  metrics  are  12%  ±  7% 
and  14%  ±  16%  for  NO  and  CO,  respectively,  with  95%  confidence. 
The  higher  confidence  interval,  i.e.  16%,  for  the  CO  case  reflects 
the  strong  variations  observed  for  such  chemical  species  during 
the  experimental  campaign.  Such  level  of  predictivity  is  extremely 
satisfying,  especially  if  we  take  into  account  the  absolute  values 
of  the  monitored  output  variables:  CO  and  NO  ranges  are  between 
6-12  ppmv  and  46-51  ppmv,  respectively.  This  means  that 
absolute  errors  in  the  predictions  of  a  few  ppm  may  lead  to  large 
errors  when  normalized  to  the  measured  value.  The  level  of 
agreement  provided  by  the  KEE-58  mechanism  is  less  satisfactory 
for  the  PHe  =  120  bar,  especially  for  the  NO  emissions,  which  are 
under  predicted  over  the  entire  range  of  the  H2  mass  fractions. 
The  average  error  metrics  (Table  5)  for  NO  and  CO  are  58%  ±  6% 
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Fig.  11.  Relative  and  global  validation  metrics  for  NO  and  CO  measurements  (dry  basis)  at  different  burner  loads,  obtained  with  the  JL  and  KEE-58  kinetic  mechanisms.  Runs 
1-6,  Table  1. 


Table  5 

Comparison  between  predicted  and  experimental  NO  and  CO  emissions  (dry  basis). 


CO 

NO 

Model  (ppmv) 

Exp  (ppmv) 

E  =  IOm 

-ye)/ye  1 

ci/% 

Model  (ppmv) 

Exp  (ppmv) 

E  =  \(ym 

-ye) /ye  1 

a  fa 

Run 

JL 

KEE 

JL 

KEE 

JL 

KEE 

JL 

KEE 

1 

284 

7.9 

12.0 

23.0 

0.33 

0.09 

2090 

66.7 

50.5 

40.4 

0.32 

0.06 

2 

280 

5.9 

6.2 

44.5 

0.04 

0.22 

2480 

44.7 

45.6 

53.4 

0.02 

0.07 

3 

423 

6.3 

5.9 

70.3 

0.06 

0.18 

2270 

46.8 

45.8 

48.6 

0.02 

0.08 

4 

268 

15.4 

17.4 

14.4 

0.12 

0.06 

1260 

25.2 

49.0 

24.7 

0.49 

0.05 

5 

291 

11.2 

7.9 

35.7 

0.41 

0.14 

1510 

15.0 

40.8 

36.0 

0.63 

0.07 

6 

300 

11.9 

6.6 

44.6 

0.81 

0.21 

1520 

14.7 

39.0 

37.9 

0.62 

0.07 

\E  /y  e\avg 

\El  /y  e\avg 

\E/ye \avg 

\Ed  /y  e\avg 

PHe  (bar) 

JL 

KEE 

PHe  (bar) 

JL 

KEE 

90 

45.9 

0.14 

0.16 

90 

47.5 

0.12 

0.07 

120 

31.6 

0.45 

0.14 

120 

32.9 

0.58 

0.06 

and  45%  ±14%  respectively,  thus  indicating  a  significant  departure 
from  the  experimental  data.  Indeed,  results  provided  by  KEE-58 
show  a  strong  improvement  with  respect  to  the  global  chemistry 
approach,  which  provides  CO  and  NO  values  higher  than  the  mea¬ 
sured  ones  by  more  than  a  order  of  magnitude  (Table  5).  The  poor 
performance  of  KEE-58  for  the  PHe  =120  bar  case  with  respect  to 
the  PHe  =  90  bar  is  probably  related  to  the  strong  approximation 
made  for  estimating  the  air  inlet  temperature  (Section  4.2).  The 
preheating  efficiency  can  be  significantly  affected  by  the  operating 
conditions,  i.e.  flow  rates,  and  considering  an  average  value  of  0.5 
can  be  inappropriate  in  some  cases,  especially  at  the  higher  loads, 
when  the  pre-heater  is  expected  to  perform  better.  Therefore,  it 
would  be  desirable  to  get  directly  such  temperature  from  the 
experimental  measurements,  to  avoid  the  consequences  of  such 
rough  approximation.  To  this  purpose,  three  additional  thermo¬ 
couples  have  been  placed  near  the  air  inlet  section  for  the  future 
experimental  campaigns. 


6.  Conclusions 

The  present  paper  presents  a  numerical  and  experimental 
investigation  of  a  micro-CHP  unit,  equipped  with  a  flameless  bur¬ 
ner  which  supplies  the  heat  required  by  a  Stirling  engine  for  the 
production  of  electrical  power.  The  activity  has  been  focused  on 
the  effect  of  adding  hydrogen  to  the  fuel  on  the  combustion  regime 
features  and,  thus,  to  assess  the  feasibility  of  flameless  combustion 
with  hydrogen-enriched  fuels. 


The  numerical  model  of  the  combustion  chamber  helped  under¬ 
standing  observed  results  and  gaining  considerable  insight  into  the 
combustion  features;  however  its  development  required  signifi¬ 
cant  efforts. 

Firstly,  geometrical  features  of  the  combustion  chamber  lead  to 
a  quite  heavy  computational  grid,  although  some  details  were  ne¬ 
glected  by  using  appropriate  procedures  (e.g.  the  heater  with 
finned  surface  was  simulated  as  an  heat  sink).  Secondly,  proper 
boundary  conditions  had  to  be  defined  in  order  to  correctly  couple 
the  oxidation  process  to  the  overall  operation  of  the  CHP  unit.  In 
particular,  the  knowledge  of  the  heat  transferred  to  the  working 
fluid  of  the  Stirling  cycle  was  required  to  properly  reproduce  the 
energy  balance  of  the  system.  To  this  purpose,  a  user-defined  sub¬ 
routine  was  developed  to  evaluate  such  energy  flux,  based  on 
information  provided  by  the  experimental  campaign,  i.e.  burner 
load,  and  by  the  numerical  simulation  itself,  i.e.  exhaust  tempera¬ 
ture  and  air  inlet  temperature.  The  external  subroutine  proved  to 
perform  satisfactory,  providing  values  of  the  heat  source  for  the 
Stirling  cycle  which  lead  to  expected  thermal  efficiencies  in  agree¬ 
ment  with  those  indicated  in  the  literature  for  similar  systems. 

Finally,  rather  detailed  kinetic  schemes  were  needed  to  capture 
flameless  combustion,  global  schemes  leading  to  unrealistic  tem¬ 
perature  and  emission  levels. 

Regarding  the  main  features  of  the  combustion  process,  results 
showed  more  uniform  temperature  distribution  in  the  combustion 
chamber  for  the  higher  loads,  being  the  system  designed  to  opti¬ 
mize  the  mixing  process  for  such  operating  conditions.  The  addi¬ 
tion  of  hydrogen  to  the  fuel  was  found  to  deeply  affect  the  flame 
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structure,  resulting  in  a  reaction  region  progressively  shifted  to¬ 
wards  the  burner  for  increasing  H2  values.  Moreover,  the  temper¬ 
ature  distribution  with  hydrogen-enriched  fuels  evidenced  the 
existence  of  a  hot  core,  characterized  by  steep  gradients,  indicating 
operation  in  an  intermediate  regime  between  flameless  oxidation 
and  a  conventional  diffusion  flame. 

The  comparison  between  the  predicted  and  measured  NO  emis¬ 
sions  indicated  a  very  good  quantitative  agreement  for  the  med¬ 
ium  burner  loads,  whereas  NO  emissions  at  full  load  were  under 
predicted,  although  the  experimental  trend  was  well  captured. 
The  reason  for  such  inaccuracies  are  mainly  to  be  found  in  the  large 
uncertainty  related  to  the  estimation  of  the  air  inlet  temperature, 
which  directly  affects  the  temperature  increase  due  to  the  oxida¬ 
tion  process  and  ultimately  influences  the  energy  balance  per¬ 
formed  by  the  user-defined  subroutine.  Following  such 
considerations,  a  new  experimental  campaign  has  been  planned 
after  a  modification  of  the  experimental  apparatus,  in  order  to 
measure  directly  the  air  temperature  after  preheating. 
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